Disclosing hidden information in the quantum Zeno effect: 
Pulsed measurement of the quantum time of arrival 
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Repeated measurements of a quantum particle to check its presence in a region of space was 
proposed long ago (G. R. Allcock, Ann. Phys. 53, 286 (1969)) as a natural way to determine the 
distribution of times of arrival at the orthogonal subspace, but the method was discarded because of 
the quantum Zeno effect: in the limit of very frequent measurements the wave function is reflected 
and remains in the original subspace. We show that by normalizing the small bits of arriving 
(removed) norm, an ideal time distribution emerges in correspondence with a classical local-kinetic- 
energy distribution. 

PACS numbers: 03.65.Xp,03.65.Ta,06.30.Ft 



I ie ha mm hi m 

A major challenge is 



I. INTRODUCTION 

The theoretical treatment of time observables is one 
of the important loose ends of quantum theory. Among 
them the time of arrival (TO A) has received much atten- 
tion lately [1 S i S, S, @, [ 
[Hj], for earlier reviews see [Tj 
to find the connection between ideal TO A distributions, 
defined for the particle in isolation, formally indepen- 
dent of the measurement method, and operational ones, 
explicitly dependent on specific measurement models and 
procedures. It is important to know, for example, what 
exactly a given operational procedure is measuring, or 
if and how a given ideal quantity may or may not be 
obtained with a particular experiment. 

Modern research on the quantum TOA owes much to 
the seminal work by Allcock Looking for an ideal 

quantum arrival-time concept, he considered that arrival- 
time measuring devices should rapidly transfer any prob- 
ability that appears at x > (x — being the "arrival 
position" ) from the incident channel into various orthogo- 
nal and time-labeled measurement channels. As a simple 
model to realize this basic feature he proposed a pulsed, 
periodic removal, at time intervals St, of the wave func- 
tion for x > 0, while the x < region would not be 
affected, see Fig. [TJ A similar particle removal would 
provide the distribution of first arrivals for an ensemble 
of classical, freely- moving particles as 5t — » 0. 

The difficulties to solve the corresponding mathemat- 
ical problem lead Allcock to study instead a different, 
continuous model with an absorbing imaginary potential 
in the right half-line, — iVoQ(x), Vq > 0, to simulate the 
detection. He argued that the two models should lead to 
similar conclusions with a time resolution of the order St 
in the chopping model or h/2Vo in the complex potential 
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FIG. 1: Schematic representation of the time of arrival mea- 
surement by periodic projection of the wave function onto the 
subspace x < at times tj, j = .., —1, 0, 1, .. separated by St. 
Figures (a) and (b) represent two instants immediately before 
and after the elimination of norm at x > 0. 



model. He then solved the Schrodinger equation with 
the complex potential and noticed that for Vo 3> |-E| ma:j; , 
where E max is a maximum relevant energy in the wave 
packet, the apparatus response vanishes, —8N/8t — > 0, 
with N = (if)\ip), because of quantum mechanical reflec- 
tion. This is one of the first discussions of the quantum 
Zeno effect, although it was not known by this name at 
that time. Eight years later Misra and Sudarshan [2(| 
generalized this result studying the passage of a system 
from one predetermined subspace to its orthogonal sub- 
space: the periodic projection method in the limit <5t — > 
was presented as a natural way to model a continuous 
measurement, however it did not lead to a time distri- 
bution of the passage but to its suppression 20]. This 
lack of a "trustworthy algorithm" to compute TOA and 
related distributions prompted them to put in doubt the 
completeness of quantum theory and has been much de- 
bated since then, see [2l|, [22|, [23[ for review. For the 
specific case of a projection onto a region of space, as in 
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the TOA procedure proposed by Allcock, several works 
have later emphasized that, in the limit of infinitely fre- 
quent measurements, the resulting "Zeno dynamics" in 
the original subspace corresponds to hard wall (Dirichlet) 
boundary conditions [2J, |25|, |26[ . 

Allcock discarded the short St (or high Vq) limit as 
useless, and examined the other limit where the mea- 
surement may be expected to be Vb-independent, Vo <C 
\E m i n \, later on. We shall not treat this limit hereafter 
but, for completeness, let us briefly recall that he got the 
current density J{t) by deconvolution of the absorbed 
norm with the apparatus response function, 



J(t) 



— (^ f (t)\(k5{x) 



■6(x)kM f (t)), (1) 



where k and x are the wavenumber and position opera- 
tors, m is the mass, and the average is computed with 
the freely moving wave function ipf. This result makes 
sense classically, but in quantum mechanics J(t) is not 
positive semidefinite even for states composed only by 
positive momenta [l7|, H||. He also got, within some ap- 
proximation, a positive distribution, n^-(t), which has 
been later known as "Kijowski's distribution" 0, [13, H3, 



n K {t) = -^ f (t)\k^5(x)k 1 ^ f {t)). 



(2) 



In this paper, instead of disregarding the Zeno limit be- 
cause of reflection, as it is customary, we shall look at 
the small amount of norm SN detected, i.e., eliminated 
in the projection process, and normalize, 



-SN/St 
hm — -. 

<5t^o 1 - N(oo) 



(3) 



In this manner, a rather simple ideal quantity will 
emerge: there is, in other words, interesting physical in- 
formation hidden behind the Zeno effect, which can be 
disclosed by normalization. To fulfill this program we 
shall put the parallelism hinted by Allcock between the 
pulsed measurement and the continuous measurement on 
a firmer, more quantitative basis. 



II. ZENO TIME DISTRIBUTION 

We shall now define formally the pulsed and continuous 
measurement models mentioned above and also an inter- 
mediate auxiliary model El that will be a useful bridge 
between the two. Ref. [28[ is followed initially although 
the analysis and conclusions will be quite different. 

The "chopping process" amounts to a periodic projec- 
tion of the wave function onto the x < region at instants 
separated by a time interval St. There is nothing here be- 
yond standard measurement theory (29| . Each chopping 
step eliminates interference terms in the density operator 
between right and left components, and the right compo- 
nent is separated from the ensamble (detected) so that it 



cannot come back to the left. The wave functions imme- 
diately after and before the projection at the instant tj, 
are related by ip(x,tj + ) — ip(x, tj_)Q(— x). 

The wave at x > may also be canceled with a 
"kicked" imaginary potential Vk — V5t Fst(t), where 
the subscript "k" stands for "kicked" and Fst(t) = 



provided 



V = -Wi = -iV @(x), 



v st > n. 



(4) 



(5) 



The general (and exact) evolution operator is obtained 
by repetition of the basic unit 



U k (0,6t) = e -iH 6t/h e -iV5t/h^ 



(6) 



where H = -{h 2 /2m)d 2 / dx 2 . 

For the continuous model, the evolution under the 
imaginary potential Q is given by 

U(0, St) = e-iiHo+VjSt/h _ e -iH St/H e -iV6t/H 

+ O(St 2 [V,H ]/h 2 ). (7) 

Comparing with Eq. © we see that the kicked and con- 
tinuous models agree when 



6t 2 \([V,Ho})\/n 2 ^1. 



(8) 



At first sight a large VoSt/h, see Eq. ([5]), seems to be 
incompatible with this condition so that the three mod- 
els would not agree [2"§| , In fact the numerical calcula- 
tions give a different result and show a better and better 
agreement between the continuous and pulsed models as 
Vq — » oo when St and Vq are linked by some predeter- 
mined (large) constant a, St — aU/Vo- Figures [2] and [3] 
illustrate this agreement: in Fig. [2] (inset) the average 
absorption time, {t) = / °° {—dN/dt)tdt/[l — A^(oo)], is 
shown versus St and IOH/Vq for a Gaussian wave packet 
sent from the left towards the origin. The main figure 
shows the time averages versus St (chopping) and fr/2V~o 
(continuous). The lines bend at high coupling because 
of reflection. The normalized absorption distribution as 
Vo — » oo was derived in [34| by making explicit use of the 
known stationary scattering wave- functions, 
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(iP f {t)\k5{x)k\^ f {t)), 



(9) 



where k^Ti is the initial average momentum. Fig. [3]shows 
for a more challenging state, a combination of two Gaus- 
sians, that this ideal distribution becomes indistinguable 
from the normalized chopping distribution when St is 
small enough. Even the minor details are reproduced, 
and differ from J and 11^, also shown. 

To understand the compatibility "miracle" of the 
inequalities ([5]) and ([5]), we apply the Robertson- 
Schrodinger (generalized uncertainty principle) relation, 



([V,H ]) 



< 2|AV7|AiJ , 



(10) 
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FIG. 2: Average absorption times evaluated from —dN/dt 
(normalized) for the projection method and the continuous 
(complex potential) model. The initial state is a minimum 
uncertainty product Gaussian for 23 Na atoms centered at 
xo = —500pm with Ax = 23.5/im and average velocity 0.365 
cm. In all numerical examples negative momentum compo- 
nents of the initial state are negligible. 
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FIG. 3: Time-of-arrival distributions: Flux J (dashed line), 
II a- (solid line), Hze.no (big sparse dots), and U ch o PP ing (for 
pulses separated by St — 0.266 ns, dotted line). The ini- 
tial wave packet is a combination ip = 2~ 1 ^ 2 (tpi + 1P2) of 
two Gaussian states for the center-of-mass motion of a single 
Caesium atom that become separately minimal uncertainty 
packets (with Axi = Ax2 = 0.021 /im, and average velocities 
(v)\ — 18.96 cm/s, (11)2 = 5.42 cm/s) at x = and t = /js. 

where A denotes the standard deviation. Since | AV7 is 
rigorously bounded at all times by Vb/2 [38| . imposing 
StVo — ah with a » 1, a sufficient condition for dynam- 
ical agreement among the models is 

V > AH . (11) 

For large Vq the packet is basically reflected by the wall 
so that AHq tends to retain its initial value and Eq. (TTTj) 
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FIG. 4: Diminishing ratio AHq/Vq with increasing Vb versus 
time for the same initial wave function as in Fig. 2 ( 23 Na 
atoms.) Eo is the initial kinetic energy. 

will be satisfied during the whole propagation. Fig. 0] 
shows the ratio AHq/Vq for three values of Vb as a func- 
tion of time. 

This implies in summary that Tlzeno = rijv, Eq. ([5]), 
a very remarkable result, which illustrates that an active 
intervention on the system dynamics may after all pro- 
vide an ideal quantity defined for the system in isolation. 



III. AN APPROXIMATE RELATION BETWEEN 
PULSED AND CONTINUOUS MEASUREMENTS 

So far we have discussed the limits St, Ti/Vq — * in or- 
der to find the corresponding time distribution. In fact a 
very simple argument relates the pulsed and continuous 
measurements approximately for finite, non-zero values of 
St and h/Vo, when they are sufficiently large to make re- 
flection negligible: the average detection time is delayed 
with respect to the ideal limit corresponding to Tlzeno as 

(*> « (t)zeno + St/2 « (t) Zeno + -£r, (12) 

ZVq 

see Fig. [2j since, once a particle is in a; > 0, -^r and 
St/ 2 are precisely the average life times in the continuous 
and discrete measuring models, respectively [39(. This 
suggests an approximate agreement between projection 
and continuous dynamics provided that the relation St « 
T1/V0 is satisfied. For large Vq, this is asymptotically not 
in contradiction with the requirement of a large a since 
Vq^ 1 — (aVb) -1 — > as Vb — > 00; in any case quantum 
reflection breaks down the linear dependence, see Fig. [5] 
A similar relation between pulsed and continuous mea- 
surements was described by Schulman (30j and has been 
tested experimentally [3l|. The simplest model in [3(| 
may be reinterpreted as a two-level atom in a resonant 
laser field, with the excited state decaying away from 
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the 2-level subspace at a rate 7 [H|, H = § (£_" 7 J- 
The relation between pulsed and continuous measure- 
ments follows by comparing the exponential decay for 
the effective 2-level Hamiltonian with Rabi frequency f2 
and excited state lifetime I/7, with the decay dynamics 
when 7 = and the system is projected every St into the 
ground state. It takes the form (301 St = — for 7/fi <C 1 
(weak driving). In our TOA model we have a different 
set of parameters but a comparison is possible by tak- 
ing into account that the imaginary potential (j4|) may 
be physically interpreted as the effective interaction for 
the ground state in the weak driving regime, for a lo- 
calized resonant laser excitation with subsequent decay, 

H = H + I °?%> )• This gives fl B % = ^, 

so that St « H/Vq becomes 



St 



(13) 



different from Schulman's relation, as it may be expected 
since the pulsed evolution depends on Q in Schulman's 
model but not in our case, where it is only driven by the 
kinetic energy Hamiltonian Hq. 



IV. DISCUSSION 

The first discussions of the Zeno effect, understood 
as the hindered passage of the system between orthogo- 
nal subspaces because of frequent instantaneous measure- 
ments, emphasized its problematic status and regarded 
it as a failure to simulate or define quantum passage-time 
distributions. We have shown here that in fact there is 
a "bright side" of the effect: by normalizing the little 
bits of norm removed at each projection step, a physical 
time distribution defined for the freely moving system 
emerges. (There are other "positive" uses of the Zeno 



effect, such as reduction of decoherence in quantum com- 
puting, see e.g. [32L l33|.) In the case of the projection 
measurements to determine the TOA, this distribution is 
given by Eq. ©. This result is fundamentally different 
from the current density or from Kijowski's TOA dis- 
tribution, (operational approaches to measure them by 
fluorescence are described in P, 0| ) and corresponds to a 
classical time-distribution of local kinetic energy density 
[3^ . [35( 1 , rather than a classical TOA distribution [Z(| • 

Experimental realizations of repeated measurements 
will rely on projections with a finite frequency and pulse 
duration that provide approximations to the ideal re- 
sult. Feasible schemes may be based on pulsed localized 
resonant-laser excitation [31J , or sweeping with a detuned 
laser (36[. A challenging open question for further re- 
search is the possibility to devise alternative (non peri- 
odic) pulse sequences to enhance the robustness of the 
Zeno effect in a TOA measurement, similar in spirit to 
the ones proposed in the context of quantum information 
processing [37J. Nevertheless, the ordinary (periodic) se- 
quence has successfully been applied in quantum optical 
realizations of the quantum Zeno effect [31[ • 

The proposed normalization method may be applied 
to other measurements as well, i.e. not only for a TOA 
of freely moving particles, but in general to first passages 
between orthogonal subspaces, and it will be interest- 
ing to find out in each case the ideal time distribution 
brought out by normalization. 
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